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ABSTRACT 

We investigate the figure rotation of dark matter halos identified in ACDM simulations. We 
find that when strict criteria are used to select suitable halos for study, 5 of the 222 halos 
identified in our z=0 simulation output undergo coherent figure rotation over a 5 hr 1 Gyr 
period. We discuss the effects of varying the selection criteria and find that pattern speeds 
for a much larger fraction of the halos can be measured when the criteria are relaxed. Pattern 
speeds measured over a 1 fe" 1 Gyr period follow a log-normal distribution, centred at fi p = 
0.25 h rad Gyr -1 with a maximum value of 0.94 h rad Gyr -1 . Over a 5 h^ 1 Gyr period, the 
average pattern speed of a halo is about 0. 1 h rad Gyr -1 and the largest pattern speed found is 
0.24 h rad Gyr . Less than half of the selected halos showed alignment between their figure 
rotation axis and minor axis, the exact fraction being somewhat dependent on how one defines 
a halo. While the pattern speeds observed are lower than those generally thought capable of 
causing spiral structure, we note that coherent figure rotation is found over very long periods 
and argue that further simulations would be required before strong conclusions about spiral 
structure in all galaxies could be drawn. We find no correlation between halo properties such 
as total mass and the pattern speed. 

Key words: galaxies: kinematics and dynamics - galaxies: halos - methods: N-body simula- 
tions - dark matter 



1 INTRODUCTION 

Initially, calculations of stellar orbits in galaxies ignored the dark 
matter halo completely. Then, as ob servational evidence for dark 
matter halos increased dGunnl fl980l) . orbit calculations involved 
modelling the halos as stationary spherical objects. Further inves- 
tigation into dark matter halos indicated that many halos we re in 
fact triaxial dDubinski & Carlberdll99ll : IWarren et al. Ill992h and 
not spherical. This w as taken into consideration in, for example, 
ISchwarzschiicidl993l) but orbits were still modelled within station- 
ary dark matter halos and th ese calculations fail to take into account 
the effect of figure rotation. Bekki & Freeman ( 2002) suggest that 
figure rotation of the dark halo would play an important role in the 
formation and evolution of the embedded galaxies. In particular, 
they suggest that figure rotation may influence the formation of 
stellar bars as well as spiral arms and warps. They also suggest 
that figure rotation could trigger star-bursts in galaxies, even at 
high redshifts. Since orbits passing near the centre of galaxies 
are affected, figure rotation c ould also have conseque nces for 
black hole growth in galaxies dde Zeeuw & Franx|ll99lh and the 
so-called 'cusp/core problem' dde Blok et alj|2001f) . Observations 
of such effects could be used to constrain current galaxy formation 
theories. 



Initial work on the subject of figure rotation was done by 
Dubinski in 1992. He found that halos in simulations which in- 
cluded tidal effects experienced a noticeable rotation. The halos 
were found to have pattern speeds ranging from 0.1 to 1.6 radians 
per Gyr. (These units are essentially equivalent to km s _1 kpc~ ). 
The small sample size and somewhat artificial implementation of 
the tidal field, meant further studies in this area were necessary. 
iPfitzneil (1999), noticed that a significant number of triaxial halos 
produced in cold dark matter simulations tended to rotate steadily 
around their minor axes. One such halo was extracted from the sim- 
ulation at late times and was allowed to evolve in isolation over 
a 5 Gyr period. During this time the halo was found to rotate as 
a solid body with con stant speed of about 1.1 radians per Gyr. 
In 1999. iBureau et alj proposed that slow figure rotation of a sur- 
rounding triaxial dark halo may be responsible for creating the spi- 
ral structure seen in NGC 2915 and encouraged further studies of 
gas di sks contained within the potentials of rotating triaxial dark 
halos. Bekk i & Freeman! J2002t) continued this work, using numer- 
ical simulations. They showed that a 'slowly' rotating triaxial dark 
halo could cause large spiral arms in an extended gas disk, and that 
the resulti ng structure is strongly d ependent on the pattern speed of 
the halo. lMasset & Bureau! (2003) furt her explored the effe cts of a 
rotating triaxial dark halo proposed by [Bureau et alj dl999h . using 
hydrodynamical simulations and concluded that, while a rotating 
triaxial halo could be responsible for causing spiral structure, the 
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required pattern speed was prohibitively large - between 5.5 and 
6.5 radians per Gyr. 

More recently, Bailin and Steinmetz (2004, hereafter BS04) 
compared the orientation of the major axis determined from the 
central sphere (radius of 0.6 x virial radius) of 317 halos selected 
from a high resolution simulation over five time intervals (a total 
period of just over 1 Gyr). They found that 278 (88 per cent) of their 
halos did indeed rotate smoothly, with pattern speeds of around 
0.15 h (where Hq = 100 h km s~ Mpc" 1 ) radians per Gyr. They 
noted that the pattern speeds observed in the inner regions of the 
undisturbed halos in their samp le were not suffici ent to cause the 
spiral patterns as discussed bv lBureau et al. 1 19991). 

We note that the simulations of Mas set & Bureaul J2003h arc 
2-D and the code used was usually applied to d isks around stars. 
When one considers the 3-D simulations of Bek ki & Freeman! 
(2002), in their figure 2f, some spiral structure is evident af- 
ter 1 G yr for much smal l er pat tern speeds than those consid- 
ered by Masset & Bureaul (2003). The slowest figure rotation in 
iBekki & Freeman! ( 120021) " is 0.77 kms kpc , which is signifi- 
cantly larger than typical values found in BS04, but it is unclear 
what would be expected in galaxies with slower pattern speeds over 
much longer periods of time. Here we extend the analysis of BS04 
to investigate whether one would expect coherent figure rotation 
over periods much longer than the 1 Gyr period they used. 

We have simulated the evolution of structure in a (50 h~ 
Mpc) 3 region of space, in a ACDM c osmology, using the N-body 
code GADGET (Springel et alj|200ll) to evolve the positions and 
velocities of 256 3 collisionless particles. Halos are identified from 
the simulation using a friends-of-friends algorithm. To explore fig- 
ure rotation, we have identified halos from simulation outputs at a 
redshift of zero and then traced them back through time. We then 
determine the principal axes of each halo using the inertia tensor. 
By following the motion of the principal axes of each halo over 
several time steps, we are able to measure the pattern speed of the 
halo. Our approach is similar to that of BS04 but we extend their 
analysis to test some of their assumptions. 

The paper is laid out as follows. In section 2, the simulations 
and analysis procedure are described. In section 3, we present the 
results of the analysis: the mass function, the spin parameter distri- 
bution, the outcome of halo tracing and the accretion constraints 
applied, the outcome of substructure elimination, the environment 
of halos selected for study, the figure rotation measurements, the 
alignment between rotation axis and the minor axis and, finally, 
the correlation between the pattern speed and halo properties such 
as halo mass. In section 4, we discuss the results. 



Parameter 


Parameter Value 


^particles 


256 3 


Box Length 


50 ft -1 Mpc 


^particle 


6.2 X 10 s ft -1 M 


Force Softening 


30 ft- 1 kpc 


Initial Redshift 


55 



Table 1. Simulation parameters 



clusters of galaxies. To set initia l particle positions and velocities, 
we us ed the COSMICS package (van de Weygae rt & Bertschinger I 
1996) with fl m =0.3, n A =0.7 and <r 8 =0.9. 



2.2 Simulations 

We ran the publicly available parallel version of GADGET on a 
cluster of eight 2 GByte machines. Doing so allowed us to follow 
the evolution of 256 3 particles. We have chosen to simulate a cube 
of length 50 h~ l Mpc. T his is sufficiently l arge to produce reli- 
able periodic simulations dPower et al.ll2003l) . The choice of par- 
ticle number and volume mean that in this simulation each dark 
m atter particle has a mass o f ~ 6x 1O 8 M0. Guided by studies 
of I Athanassoula et al. I < f2OO0l) we have used a force softening of 30 
/i _1 kpc to prevent numerical instability. The parameters chosen 
for our simulation are summarised in TableQ] 



2.3 Analysis procedure 

We have identified halos (groups of ten or more particles) in the 
redshift zero snapshot using FoF, a simple friends-of-friends group 
finder for N-body simulations developed by the simulation group 
at the University of Washington. We have used the common con- 
vention of setting our linking length to be 0.2 times the mean inter- 
particle spacing. We obtained results for two different definitions of 
a halo. Firstly, we used all particles identified by the group finder 
and relied on our technique for substructure elimination to exclude 
groups which contained 'sub-halos'. We also tried the halo defini- 
tion of BS04, which selects a central sphere of particles within a 
group. Once we have defined which particles are in a halo, we use 
the positions of these particles to calculate the inertia tensor, and 
hence the principal axes, of the halo of interest at each time step. 



2 METHODOLOGY 



2.1 Initial conditions 



We assume the standard model for structure formation: that density 
fluctuations are gaussian and the initial power spectrum is given 
by a power law, Po oc k n , where n ~ 1. As the universe moves 
from a radiation-dominated to a matter-dominated phase and 
fluctuations begin to grow, the power spectrum bec omes modified 
in a w ay that depends on the density of dark matter. Barde enet al.l 
d 19861) give the transfer function T (where P(k) = Po(k) T) for a 
universe that is dominated by cold dark matter. The normalisation 
of the spectrum can be obtained from observations of the Cosmic 
Microwave Background or from the measured abundance of 



2.3.1 Calculating the inertia tensor 

For figure rotation analysis, we model each halo as a rigid body, 
and calculate the inertia tensor using 



hj = m a 



k 



The eigenvalues and eigenvectors of this tensor are then extracted 
using Jacobi Transformations. We have also performed the figure 
rotatio n analysis using the modified inertia tensor suggested by 
BS04, which artificially weights central particles in halos. 
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2.3.2 Tracing halos through time and accretion constraints 

To follow halo properties with redshift, we need to associate each 
halo with its progenitor halos at a previous timestep. This is 
done w ith merger trees using code developed by Kauffmann et al. 

The algorithm begins at the first output time, Zi, at which 
a halo (of 10 or more particles) can be found. It then steps through 
all halos identified at the next output time, Zi+i, searching for pro- 
genitor halos in the z% output. A halo at redshift Zi is said to be a 
progenitor of a halo at redshift Zi+\ if, and only if, two conditions 
are satisfied: 

(i) A progenitor halo at z; must contribute at least half of its 
particles to the halo at Zi+\ and 

(ii) the central (most gravitationally bound) particle of the pro- 
genitor halo at Zi must be contained within the halo at Zj+i 

We then step through halos in the Zi+2 output, searching for 
progenitors in the Zi+i timestep. This continues until the output 
time of z=0 is reached. 

Once a merger tree has been established, we exclude halos 
which have accreted a significant fraction of their mass over the 
time period considered. In the "halo matching" technique of BS04, 
they excluded halos which accreted more than 15 per cent of their 
mass over a 1 ft -1 Gyr period; that is, to be considered for figure 
rotation analysis, a halo at z=0 must contain at least 85 per cent of 
the particles contained by the largest progenitor halo identified in 
the timestep which is 1 h~ x Gyr earlier. 

Since we extend our analysis to consider halos over a much 
longer period of time, we have explored the effects of relaxing this 
limit over the 3 and 5 hT 1 Gyr period and include halos which ac- 
crete as much as 25 per cent of their mass over 3ft" 1 Gyr and 30 
per cent over 5 hT 1 Gyr. We then consider the consequences of re- 
laxing the constraint even further in our 5 hT 1 Gyr analysis: instead 
of comparing particles in a halo at z=0 with the largest progenitor 
halo identified 5 h^ 1 Gyr before, we impose a constraint on pro- 
genitor halos at each successive timestep. Thus, at each timestep, 
halos must contain at least 70 per cent of the particles contained by 
the largest progenitor halo in the previous timestep. 

2.3.3 Measuring the figure rotation of a halo 

To determine the figure rotation of the halo we then followed a 
plane fitting approach suggested by BS04, where a detailed descrip- 
tion of the method is given. Briefly, the method involves fitting a 
plane to the principal axes of the halo of interest determined over 
several time steps and then solving for the plane, z = ax + by, that 
best fits these axes. Once the plane has been determined, measuring 
the amount of rotation in the plane is straight forward. The princi- 
pal axis is projected onto the plane at each time step and the angle 
between the projected principal axes over a time interval gives the 
figure rotation of the halo. We then used linear regression to find the 
best-fitting linear relation for this figure rotation, the slope giving 
the pattern speed of the halo. 

We have also measured the figure rotatio n of several halos us- 
ing the method described by Dubinski ( 1992) and found that these 
results are similar to the ones we obtain using the plane fitting ap- 
proach. 



3 RESULTS 

Many of the results from this section are summarized in Table 2. 



o 

Q-3.5 
O 

7 3 
£2.5 

o 

D) 
O 

5 2 

n 

T3 



11 11.5 12 12.5 13 13.5 14 

log10(M//T M Q ) 

Figure 1. The mass function. The histogram shows the number of halos as 
a function of the halo mass in the simulated region at redshift zero and the 
line shows the fit obtained by Jenkins et al. (2001). 



3.1 Mass function 

The mass function is a measure of the n umber of halos as a func- 
tion of their mass. Ijenkins et al] J200lh were able to predict the 
mass distribution of dark halos expected in a Cold Dark Matter uni- 
verse by combining the results from sev eral N-body simulations. 
The fit obtained by Uenkins et al ] d200lh shows good agreement 
with the theoretical pre dictions of the Sheth and Tormen formula 
( Sheth &TormerJl999l) . The mass function obtained from our sim- 
ulation is shown in FigureQ] This figure is a log plot of the number 
of halos as a function of their mass in the (50 h' 1 Mpc) 3 region 
of our simulation, at a redshift of zero . We have fo u nd th at our 
mass function agrees well with that of I Jenkins et al.l j200lh over 
the range of halos simulated. 



3.2 Angular momentum 

We use the comput ationally convenient expression for the spin pa- 
rameter, defined by Bullock et al. (200 lh as 



V2MVR 

to calculate the spin parameter of each of our simulated halos. This 
expression describes the spin parameter for a sphere of radius R 
(taken to be the virial radius), where M is the mass contained within 
the sphere, J is the angular momentum of the sphere, and V is the 
circular velocity at the radius R. Figure [2] shows the distribution of 
the spin parameter found in our simulation. The spin parameter was 
calculated for all simulated halos which have more than 200 parti- 
cles (M > 10 11 /i _1 M ).The curve is the log-normal distribution 
given by 

P(X') = exp (- ^ (2) 

with best-fitting values of X' = 0.034 ± 0.003 and o = 
0.62 ± 0.08. We have found that our best-fittin g value for An and 
a are cons istent with the value s obtained by iBarnes & Efstathio ul 
dl987l) and lBullock etai] feOOlh . 
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1 ft- 1 Gyr 


3 h- x Gyr 


5 h.- 1 Gyr 


5 h Gyr relaxed 


Halos traced 


218 


201 


191 


191 


Halos after accretion cut 


183 


62 


13 


74 


Halos after substructure cut 


115 


38 


5 


48 


Halos after BS04 substructure cut 


89 


55 


12 


70 


Coherent Figure Rotation 


61% 


71% 


100% 


73% 


fj, (fitted) 


-0.60 ±0.03 


-0.95 ±0.04 


-1.03 


-0.94 ±0.02 


fi (rad/Gyr) 


0.25h 


0.1 lh 


0.09h 


0.1 lh 


a 


0.37 ± 0.03 


0.27 ± 0.04 




0.21 ± 0.02 



Table 2. Summary of results. Values given in the first four rows correspond to the number of halos, of the original 222, which remain after various cuts are 
made. The accretion cut is applied before the substructure cut. We give results for two different methods of doing the substructure cut. Values given in the last 
three rows are for the first substructure cut, not that of BS04. The fitted value of fi is also given in units of radians per gigayear. Only the mean is given for the 
5 h~ 1 Gyr sample since it only includes 5 halos. 
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Figure 2. This figure shows the distribution of the dimensionless spin pa- 
rameter A' for all simulated halos containing at least 200 particles, together 
with the log-normal distribution given by equation [2] with fitted parameters 
\' = 0.034 ± 0.003 and a = 0.62 ± 0.08. 

3.3 Halo tracing and accretion constraint results 

BS04 used a computationally intensive bootstrapping technique to 
show that the error in the determination of a halo's principal axes 
depends on the number of particles contained within that halo and 
the intrinsic shape of the halo. They also showed that it was not 
possible to measure the principal axes of halos with the precision 
required for figure rotation measurements if the halos contained 
fewer than 4000 particles. Guided by their work, we identified the 
222 halos in our simulation at a redshift of zero, which contain at 
least 4000 particles (that is, a mass of at least 2.5 x 10 12 ft -1 M Q ). 
BS04 also showed that errors are large for halos which are almost 
oblate but this only affects a small number of halos in our case. 
We are able to trace 218, 201, 191 halos over the 1, 3, 5 ft" 1 Gyr 
periods respectively. 

After imposing constraints on the acceptable amount of ac- 
cretion over the time period considered, as discussed in !2.3.2l we 
are left with samples of 183, 62 and 13 halos over the 1, 3, 5 ft" 1 
Gyr periods respectively. When we use the relaxed 5ft" 1 Gyr con- 
straints, ie. impose a constraint at each timestep, rather than a con- 
straint over the entire period, we obtain 74 halos over 5 ft -1 Gyr. 

3.4 Eliminating halos with substructure 

We are concerned with the figure rotation of undisturbed halos 
and, therefore, needed to eliminate all halos exhibiting a signifi- 
cant amount of substructure. Identifying substructure is somewhat 
subjective and we compared various methods. The first method in- 
volved projecting the distribution of particles along three perpen- 



dicular axes and considering each axis individually. We expect the 
mass distribution of an undisturbed halo to be smooth when con- 
sidered along any axis and the presence of a secondary peak in this 
distribution thus provides evidence of substructure. We calculated 
the mass contained in bins moving out from the centre of the halo, 
where the centre is defined as the most bound particle. If any bin 
contained more mass than an inner bin we compared the mass of 
that bin to the total mass of the halo. After visually inspecting the 
distributions, and some experimentation, we chose to use 25 bins 
and eliminate halos if an outer bin contained more than 10 per cent 
of the total halo mass. Applying this to our sample of halos over a 
period of 1 ft -1 Gyr, we eliminated 68 halos on the basis of sub- 
structure. 24 and 8 halos needed to be excluded from the sample 
over the 3 and 5 ft" 1 Gyr periods respectively. For our larger 5 ft -1 
Gyr sample, where we constrain common particles between suc- 
cessive timesteps, 26 of the 74 halos need to be eliminated. 

We also tried using the method of lBS04l Their method in- 
volves excluding halos if the calculated pattern speed varies signif- 
icantly when one changes the size of the sphere used to define a 
halo. The BS04 method only considers particles within the central 
sphere of the halo and substructure outside of this will not cause 
a halo to be eliminated. Using this approach we found that 94, 7 
and 1 of the 183, 62 and 13 halos we considered would have been 
eliminated. For the relaxed constraint over 5 ft" 1 Gyr, 4 of the 74 
halos need to be eliminated. We note that our method was generally 
more conservative over longer time periods (over the 5 ft" 1 Gyr pe- 
riod we rejected 35 per cent of halos while 5 per cent would have 
been eliminated using their method) but there were a few halos 
that would have been excluded using their method and were not 
excluded using ours. 



3.5 Mass density in surrounding environment 

In order to explore the environment of our halos, we have consid- 
ered the mass density of the area surrounding the halos at z=0. To 
do this we have, for each halo, found the number of particles within 
5ft" 1 Mpc of the halo. The sum of the surrounding mass is then 
divided by the volume of a sphere with radius of 5 ft" 1 Mpc. The 
results are shown in Figure[3]where we plot the fraction of halos as 
a function of the log of their mass density. The top figure shows the 
distribution of mass density of all the halos identified from the sim- 
ulation, while the bottom histogram shows the distribution of our 
sample of undisturbed halos. While we lose a few of the most mas- 
sive halos in our selection, the undisturbed halos appear to occur 
in an environment which is sligh tly more dense than a verage. This 
is consistent with the results o f iGottlober etal] feOOlh who found 
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Figure 3. The mass density in surrounding environment. These histograms 
show the fraction of halos versus the log of the mass density of the sur- 
rounding environment. Top: The surrounding environment of all simulated 
halos. Bottom: The environment of the halos selected for 3 Gyr analysis. 

that, for z < 1, the merger rate of cluster halos to be 3 times lower 
than that of isolated halos and twice as low as halos that end up in 
groups. With this in mind, we expect that undisturbed halos would 
occur in a more dense environment. 



3.6 Figure rotation 

To determine the pattern speed of a halo, we have followed each 
halo over a period of 5 h^ 1 Gyr (or 1 or 3 h~ Gyr if the halo does 
not satisfy selection criteria over longer periods). Using linear re- 
gression, we have found the best-fitting linear relation for the figure 
rotation of the halos as a function of time. The pattern speed of the 
halo is given by the slope of the linear fit. We have calculated the 
one-sigma limit of the slope, and we have taken this to be the er- 
ror in the pattern speed. Figure [4] shows the figure rotation of one 
of our halos that has remained unaffected by substructure over a 5 
h^ 1 Gyr period. Results are similar for other halos. A linear fit is 
adequate over the whole period and there is no indication that the 
pattern speeds are systematically changing over time. 

Figure [5] shows the pattern speeds obtained for our halos. 
The 115 halos observed over 1 h' 1 Gyr are shown in the top plot, 
the middle plot shows the 38 halos observed over 3 IT 1 Gyr, the 
bottom plot shows the halos considered to be undisturbed over 5 
h^ 1 Gyr period using relaxed accretion constraints (boxed crosses 
mark the pattern speeds and errors of the 5 halos that meet the 
stricter accretion criteria similar to that used in the top and middle 
plots and in BS04). The x-axis is the estimated error in pattern 
speeds. The dashed line represents the points at which the observed 
pattern speed is equal to the estimated error. We have also shown 
the points at which the pattern speed is equal to twice the estimated 
error (dotted line). Using a 2a cutoff for halos to be considered 
rotating coherently, 63 per cent of the halos in the 1 /i -1 Gyr 




Time {h' gigayears) 

Figure 4. Pattern speed of a single halo measured over 5 h~ Gyr. Each 
point represents the angle (in radians) that the figure has rotated over the 
time interval given. The slope of this plot (obtained by linear regression) 
gives the pattern speed of the halo. 



sample, 71 per cent of the halos in the 3 h^ 1 Gyr and 100(73) per 
cent of the halos in the 5 hT 1 Gyr(relaxed) sample exhibit coherent 
figure rotation. 

We found that the distribution of pattern speeds of the undis- 
turbed halos from our simulation was fairly well fit by the 
log-normal distribution 



P(fi„) = 




Figure [6] shows the distribution of pattern speeds for the halos 
considered, together with the log-normal distribution (equation[3](. 
Best-fitting values to this curve are given in Table 2. For the undis- 
turbed halos over 5 ft -1 Gyr the fastest figure rotation detected was 
0.13 h radians per Gyr and the average pattern speed of these halos 
was 0.09 h radians per Gyr. The maximum figure rotation for the 
halos with relaxed contraints over 5 h^ 1 is 0.24 h radians per Gyr. 
Adjusting our cutoff to include only those halos with pattern speeds 
greater than three times the error, does not significantly change the 
pattern speed distribution. 



3.7 Alignment between figure rotation axis and minor axis 

iDubinskil dl992h . IPfitznerl Jl999h and lBS04l all noted that for the 
majority of their simulated halos the major axis appeared to rotate 
around the minor axis. To explore this, we have determined the dot 
product between the figure rotation axis and the minor axis of each 
of the halos in our sample. Figure [7] shows our results when we 
use a standard inertia tensor and do not limit ourselves to central 
particles in the definition of the halo. It appears that these axes are 
well-aligned for at least a third of the halos. We note, however, that 
the fraction of halos exhibiting this alignment is smaller than in 
other studies. If we consider only the spherical centre of the halos, 
the percentage of halos exhibiting this alignment increases to ~ 50 
per ce nt which is more consistent with the 58 per cent found by 
BS04. Using a modified inertia tensor also slightly increases the 
fraction of halos found to be rotating around the minor axis. 



6 S.E.Bryan and C.M.Cress 




Error in H p (h radians gigayear" ) 



log 10 £lp (h radians gigayear" 




Error in n p (h radians gigayear" ) 



log 10 il p (h radians gigayear ) 



Error in H p (h radians gigayear" ) 



m 

-1.5 -1 -0.5 

log 1 £l p [h radians gigayear" 



Figure 5. Measured pattern speeds of the simulated halos versus the esti- 
mated error in pattern speed. The dashed line represents the points at which 
the pattern speed is equal to the estimated error (la limit), the dotted line 
represents the points at which the pattern speed is equal to twice the error. 
Halos with pattern speeds less than the twice the error are not considered 
to be rotating coherently. Top: Halos considered over 1 h~ l Gyr period. 
Centre: Halos considered over 3h~ 1 Gyr. Bottom: Halos considered over 
5 h^ 1 Gyr period using relaxed accretion constraints. Boxed crosses mark 
the pattern speeds and errors of the 5 halos that meet the stricter accretion 
criteria similar to that used in the top and middle plots and in BS04). 

3.8 Effect of halo properties on pattern speed 

We have investigated the correlation of pattern speed with halo 
properties such as mass, density of the surrounding environment 
and the amount of matter accreted over 5 hT 1 Gyr (fraction of par- 
ticles common to progenitor and current halos). We plot results in 
Figure [8] demonstrating that pattern speed is not correlated with 
any of these halo properties. The final plot demonstrates that our 
results are not sensitive to the the accretion cut values we use. 



4 DISCUSSION AND CONCLUSIONS 

We have found that the mass function determined fr om our sim- 
ulation is well fit by the mass function obtained by Jenki ns et al.l 



Figure 6. Pattern speed distribution for rotating halos. This histogram 
shows the distribution of the pattern speeds calculated from our undisturbed 
halos. The line is the log-normal distribution (equation [3)- Top: Halos Ro- 
tating over 1 h Gyr (best-fitting values of fi = —0.60 ± 0.03 and rr = 
0.37 ± 0.03). Centre: Halos Rotating over 3 h' 1 Gyr (/i = -0.95 ± 0.04 
and a = 0.27 ± 0.04). Bottom: Halos Rotating over 5 h~ 1 Gyr (relaxed 
constraints) (fj, = -0.94 ± 0.02 and a = 0.21 ± 0.02). 

using N-body simulations. The distribution of the spin pa- 
rameter found from our simulation was well fit by a log-normal dis- 
tribution (equation [2} with a best-fitti ng value of \' n = 0.034. This 
agre es well with the val ues found by [Barnes & Efstathioul f 19871) 
and lBullock et alj feOOll) . Both of these tests suggest that the simu- 
lations we have run are reliable, agreeing with many previous sim- 
ulations as well as with theoretical predictions. 

We selected halos which had not accreted a large fraction of 
their material and which did not show evidence for significant sub- 
structure and found that 63 per cent of the 115 halos considered 
over 1 h~ Gyr did exhibit coherent figure rotation, where 'co- 
herent' is arbitrarily defined to be pattern speeds greater than two 
times the estimated error. The distribution of the pattern speeds 
over 1 h' 1 Gyr was well-fit by a log-normal distribution centred at 
log(fip) = -0.60 ± 0.01 (ftp = 0.24 h radians per Gyr). The pat- 
tern speeds of these halos, determined over al/i -1 Gyr period, are 
co nsistent w ith the lower end of the 0.1-1.6 radians per Gyr found 
bv lDubinskil Jl992h . but slower than the 1.1 radians per Gyr de- 
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Figure 7. Alignment between figure rotation axis and minor axis. This fig- 
ure shows the dot product of the figure rotation axis and the minor axis. 
Top: axes calculated for the 115 halos in our sample using the whole halo 
and the 'standard' inertia tensor. Middle: axes calculated for the 73 halos 
found to undergo coherent rotation. Bottom: axes (for 115 halos) calculated 
using a central sphere of 0.6 times the virial radius and a modified inertia 
tensor. 



tected bv lPfitzned fl999h . Our distribution of patter n speeds peaks 
at slightly higher values than those found by BS04, but the differ- 
ence is not very significant. It appears that limiting oneself to the 
central sphere of a group and using a modified inertia tensor do not 
significantly influence the results. The accretion cut and the sub- 
structure cut affect the number of halos considered, but the pattern 
speed distributions are not very sensitive to the details. 

The pattern speed distributions appear to shift to lower speeds 
as one looks over longer periods. The difference between the 3 
h^ 1 Gyr and 5 h^ 1 Gyr distributions is not significant but the dif- 
ference between the 1 ft -1 Gyr and 5 h^ 1 Gyr distributions is sig- 
nificant. This result holds no matter what halo definition, accre- 
tion cut or method of substructure removal is used. Presumably, 
studying the halos over a longer time period allows us to de- 
tect smaller figure rotations above the noise but, over 5 Gyr, 
there are no halos with the larger figure rotations seen over 1 
/i Gyr, suggesting that accretion and/or mergers allow for in- 
creased pattern s peeds. In general, the s peeds found are much 
slower than those Masset & Bureau! f2003h claim are required for 
spiral structure in NGC2915 but, given the long periods of time 
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Figure 8. Pattern speed versus halo properties. We have compared pattern 
speeds to (from top to bottom) the mass density in the surrounding envi- 
ronment, the halo mass and the fraction of current halo particles present 
in progenitor halo. The pattern speeds are calculated over the 5 ft. -1 Gyr 
period. The halo mass and mass density in surrounding environment are de- 
termined at z=0. Boxed crosses are halos which meet the stricter accretion 
criteria ofBS04. 



over w hich these halos show coh erent figure rotation and the re- 
sults of iBekki & Freeman] (2002), we suggest that further hydro- 
simulations would be required before strong conclusions can be 
drawn for all galaxies. 

The sample of undisturbed halos we have chosen for analysis 
are generally found to be in more dense environments than the av- 
erage halos in the sim ulation. This is consistent with the results of 
iGottlober et alj d200lh who suggest that, at low redshifts (z< 1), the 
merger rate is lower in denser environments. 

In studying the relationship between the figure rotation axis 
and the minor axis, we found alignment for many of the halos but a 
significant fraction which were not very well aligned. The two axes 
are within 25° of each other in 30 per cent of the halos, using the 
"whole-group" definition of a halo. If we restricted our anal ysis to 
the central region and used the modified inertia tensor, as did BS04, 
our results are more similar to theirs. 
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We have found no correlation between halo properties such as 
halo mass, or environment, and the pattern speed. 
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